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Abstract 

The STAR experiment at RHIC has a unique capability of measuring identified hadrons over a wide 
range of pseudorapidity (77), transverse momentum (pr), and azimuthal angle (<f>) acceptance. The 
data collected ( y/s^N = 7.7, 11.5, and 39 GeV) in its beam energy scan (BES) program provide a 
chance to investigate the final hadronic state freeze-out conditions of ultrarelativistic Au+Au col- 
lisions. The particle ratios are used to compare to a statistical model calculation using both grand 
canonical and strangeness canonical ensembles to extract the chemical freeze-out parameters. The 
Pt distributions are fitted to calculations using a blast-wave model to obtain the kinetic freeze-out 
parameters. We discuss the centrality dependence of the extracted chemical and kinetic freeze-out 
parameters at these lower energies. 



1. Introduction 

One of the early goals of heavy-ion collisions at Relativistic Heavy Ion Collider (RHIC) is to 
establish the existence of a new state of matter which is called quark-gluon plasma (QGP) [1]. Ac- 
cording to lattice quantum chromodynamics calculations (QCD), the theory of strong interactions, 
QGP occurs when a sufficiently high temperature and high energy density (« 1 GeV/fm 3 ) is reached. 
The partonic system is transformed into hadronic matter as the system cools to lower temperatures. 
The phase diagram of QCD is in general characterized by two quantities, the temperature (T) and the 
baryon chemical potential (ps) or the (net) baryon density (n#). The phase diagram should contain 
information about the phase boundary that separates the QGP and hadronic phases [2]. Lattice QCD 
finds a rapid, but smooth crossover transition from hadron gas to QGP at vanishing baryon chemical 
potential and large temperature T, while various models predict a strong, first-order phase transition 
at large pg- If this is the case, then there should be a critical point at intermediate values in the (T, pg 
) plane where the transition changes from a smooth crossover to a first order [3]. The BES program 
at RHIC is carried out using several center-of-mass energies of colliding nuclei, to explore the above 
aspects of this QCD phase diagram. 



The constituents of the hot and dense medium produced during a heavy-ion collision interact 
with each other by inelastic and elastic collisions and it evolves into a state of free particles. This 
process of hadron decoupling is called freeze-out. Two kinds of freeze-out are found: chemical 
freeze-out (T c /,) when inelastic collisions cease and the particle yields become fixed; thermal (ki- 
netic) freeze-out (T^„) when elastic collisions cease and particle transverse momenta (pj) spectra 
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Figure 1 : (a) Simultaneous blast-wave fit of charged pion, kaon and proton for 0-5% centrality in Au+Au collisions at yijviv = 
11.5 GeV. (b) Variation of T K „ with </3) for different centralities at 7.7, 11.5, 39, 62.4, and 200 GeV. The 62.4 and 

200 GeV results are taken from Ref.[4]. The errors shown here are the statistical and systematic errors added in quadrature. 



get fixed. 

We present a study of the centrality dependence of hadronic freeze-out parameters in Au+Au col- 
lisions at mid-rapidity for s/snn = 7.7, 11.5, and 39 GeV measured by the STAR experiment. To 
extract chemical freeze-out parameters we use a statistical thermal model (THERMUS) [5] where 
we fit experimental particle ratios using a grand canonical ensemble (GCE) approach with the in- 
clusion of a strangeness saturation factor (75 ) and also with strangeness canonical ensemble (SCE) 
where strangeness quantum number is conserved exactly. In this study we have used mid-rapidity 
particle ratios that include the pions (ti + , kaons (K + , K~), protons (p, p), K°s, Lambdas (A, A) 
and Cascades (E~, H + ) [6, 7]. The chemical freeze-out parameters extracted are T c /,, fig, p-s and js 
. The kinetic freeze-out parameters are determined from the blast-wave model (BW) [8] fits to the 
Pt spectra of n, K and p [6], The main kinetic freeze-out parameters extracted are Ty„ and average 
flow velocity {(fi)). 



2. Results 

2.1. Kinetic freeze-out 

At the kinetic freeze-out, elastic collisions among the particles stop and the spectral shape of the 
particles get fixed. Kinetic freeze-out parameters are obtained using blast-wave model by doing the 
simultaneous fits of n, K, and p transverse momentum spectra. The BW model describes the spectral 
shapes assuming a locally thermalized source with a common transverse flow velocity field. It has 
been successfully used to describe pj spectra with three parameters - T&„, (J3), and the exponent in 
the flow velocity profile n at 62.4 and 200 GeV [4, 9]. Figure 1(a) shows the simultaneous blast-wave 
fit of n, K, p and the corresponding antiparticles for - 5% centrality in Au+Au collisions at ysMV 
= 11.5 GeV. The variation of T kin as a function of </3> at V*aw= 7.7, 1 1 .5, 39, 62.4 and 200 GeV is 
shown in Fig 1(b). The 62.4 and 200 GeV results are taken from the Ref.[4]. The T^,, decreases from 
peripheral to central collisions. It also decreases with increasing collision energy. The (J3) increases 
with increase of energy as well as collision centrality. So, higher value of T^ corresponds to lower 
value of (J3) and vice-versa. The errors shown are the quadratic sum of statistical and systematic 
errors. 
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Figure 2: Statistical thermal model [5] fits to experimental mid-rapidity particle ratios for 0-5% centrality in Au+Au 7.7 GeV 
energy. 
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Figure 3: Results from a statistical thermal model fit [5] for Au+Au 7.7,11.5, and 39 GeV. Chemical freeze-out temperature, 
baryon chemical potential, and strangeness saturation factor are shown as a function of N pan . The 200 GeV results are taken 
from Ref. [10]. 



2.2. Chemical freeze-out 

At the chemical freeze-out, inelastic collisions among the particles stop and hadron yields get 
fixed. A statistical thermal model (THERMUS) was used to fit mid-rapidity particle ratios including 
yields of n, K, p, K°s, A and 3 measued in Au+Au collisions at v^aw= 7.7, 11.5, and 39 GeV. 
Although the particle ratios are obtained at y=0, the measurements of yields for n, K, p are for 
|y| < 0.1 and those for K° s , A , H are for |y| < 0.5. The errors on particle ratios including yields of 
7T, K, p, K°s, A, and 3, are the quadratic sum of statistical and systematic uncertainties. Pion yields 
have been corrected for feed-down from K°s weak decays. Proton yields have not been corrected 
for feed-down contributions. The A yields have been corrected for the feed-down contributions from 
H and 3° weak decays [7]. In the framework of this model, the particle yield ratios can be described 
by a set of parameters such as T e /,, ps, P-s and js ■ The errors on freeze-out parameters are obtained 
from THERMUS model. Considering grand canonical formulation of this model we have studied 
the centrality and energy dependence of the freeze-out parameters. 

Figure 2 shows the statistical model fits to experimental particle ratios for 0-5% centrality in 
Au+Au collisions at yfsmj = 7.7 GeV. The data and model matches very well, except for A/ A. Fig- 
ure 3 (a) shows the T c /, increases with increase of collision energy. Figure 3(b) shows ps decreases 
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Figure 4: Results from a statistical thermal model fit [5] for Au+Au 7.7 GeV in GCE and SCE. Chemical freeze-out temper- 
ature, baryon chemical potential and chisquare per degrees of freedom are shown as a function of N pan . 



with increasing collision energy. From peripheral to central collisions, /jb increases at lower ener- 
gies. We observe a centrality dependence of chemical freeze-out curve (T c /, vs. fig) at BES energies 
which was not observed at higher energies like Au+Au 200 GeV [10]. Figure 3 (c) shows that the 
ys increases from peripheral to central collisions for all energies. 

In contrast to GCE where all quantum numbers are conserved on an average, THERMUS model also 
allows for a strangeness canonical ensemble where only the strangeness quantum number is required 
to be conserved exactly where as baryon and charge quantum numbers are conserved on an average. 
Figure 4 (a) and (b) indicate that in peripheral collisions, T c /, and fig follow a different behavior 
in GCE and SCE at -\/Jnn = 7.7 GeV. We observe a higher in SCE in comparison to GCE at 
pheripheral collisions. Further systematic investigations are ongoing towards a more quantitative 
analysis for all the BES energies. 



3. Summary 

The new measurement at BES energies 7.7, 1 1 .5, and 39 GeV at RHIC extends the hb range from 
20 to 400 MeV in the QCD phase diagram. Kinetic freeze-out parameters are obtained using the 
measured particle spectra and a BW model. For all the beam energies studied, the central collisions 
are characterized by a lower and larger (fi) while the peripheral collisions are found to have a 
higher T^,,, and smaller (J3). Chemical freeze-out parameters are obtained using the measured particle 
ratios and a THERMUS model. We have observed a centrality dependence of the chemical freeze- 
out parameters at the lower energies. We have observed different behavior of chemical freeze-out 
parameters (T c h> Hb) for peripheral collisions in GCE and SCE. 



References 
References 

[1] J. Adams et al. (STAR Collaborati on), Nucl . Phys. A757, 102 (2005). 

[2] P Braun-Munzinger et al. larXiv:1101.3167l 2011; B.Mohanty, Nucl. Phys. A 830, 899C (2009). 

[3] S. Gupta et al., Science 332, 1525 (201 1); E. S. Bowman and J. I. Kapusta, Phys. Rev. C 79, 015202 (2009). 

[4] B. I. Abelev et al. (STAR Collaboration), Phys. Rev. C 79, 034909 (2009). 

[5] J. Cleymans et al., Computer Physics Communications, 180, 84 (2009). 

[6] L. Kumar (STAR collaboration). larXiv:12014203l (2012). J. Phys. G: Nucl. Part. Phys. 38, 124145 (201 1). 
[7] X. Zhu (STAR Collaboration), Acta Phys. Polon. B Proc. Supp. 5 (2012) 213-218. 
[8] E. Schnedermann et al., Phys. Rev. C 48, 2462 (1993). 

4 



[9] B. I. Abelev et al. (STAR Collaboration), Phys. Rev. C 81, 02491 1 (2010). 
[10] M. M. Aggarwal et al. (STAR Collaboration), Phys. Rev. C 83, 024901 (2011). 



5 



